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• ■ 1 Introduction 

O 

CN I The theory of multiple orthogonal polynomials, rooted in Hermite-Pade approximation schemes 

[U [2], is currently undergoing development O HI O [6]. These are polynomials of a single 

variable, labeled by a multi-index and orthogonal for a number of measures. Three of us have 

^ \ recently shown [7] that the multiple Charlier polynomials arise in the construction of the common 

eigenstates of a set of non-Hermitian oscillator Hamiltonians. In keeping with the significant 

connection between algebras and special functions, this provided a model from which many 

properties of the multiple Charlier polynomials could be derived and interpreted while giving 

a physical context for them. The purpose of this paper is to similarly describe another family 

of multiple orthogonal polynomials, namely that of the multiple Meixner polynomials of the 

first kind. As shall be seen, these are related to the Lie algebra of diffeomorphisms in C 

whose generators are realized in terms of the creation and annihilation oscillators of independent 

harmonic oscillators. This will allow to obtain their step relations and their generating function. 

2 Multiple Meixner polynomials of the first kind 

The monic r-multiple Meixner polynomials of the first kind M"'^{x) are labeled by a multi-index 
n = (ni, • • • , n^) E N^ with length \n\ = ni + ■ ■ ■ + Ur- They are orthogonal with respect to the 
r discrete measures 

oo ^ 

^. = V(/3)..^5., i = l,---,r, (2.1) 

x=0 
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where /3 > 0, < q < 1 for all i and Cj ^ Cj whenever i ^ j; that is, they satisfy 

OO ,j. 

5^m|'^>V(/3).^ = 0, (2.2) 



x=0 

for j = 0, ••• ,nj — l,i = I,--- ,r. As usual, the Pochhammer symbol stands for {(3)x = 
j3{f3 + 1) • • • (/3 + X — 1). The parameters ci, • • • , c^ are collectively represented by c. In the 
following, we shall often omit the suffices /? and c whenever the context is unambiguous. 

Let ei = (0, • • • , 0, 1, 0, • • • ,0) denote the i-th standard unit vector in r dimensions with 1 in 
the z-th entry. The multiple Meixner polynomials of the first kind obey [5] r recurrence relations 
of nearest-neighbor form that read as follows: 



xMfi{x) = Mfi+g^{x) + 



Ci /o , i-i\ , V^ "-fc 



fe=l 



MJx) 



'-n\ 



(2.3) 
+ X] n -c )2 "fc(^+ 1^1 ~ l)^n-efc(a;), i = !,••• ,r. 



fc=i 

It should be observed, that if one subtracts the above relations pair-wise, one finds that the 
polynomials M^(x) verify as a consequence of (j2.3p : 

M„^+ei(x) - Mn+e, = (/3 + 1^1) (Y3^y^rr^^"(^) (2.4) 

for all i,j G {1, • • • ,r}. 

3 Harmonic oscillator operators and states 

The construction of the algebraic model for the polynomials M^' will involve the operators and 
eigenstates of r independent quantum harmonic oscillators . We here briefly review basic facts 
in this connection to establish notation. 

The annihilation and creation operators aj, af (resp.), i = 1, • • • ,r, satisfy the commutation 
relations 

[ai,aj] = [a+,a+] = 0, [ai,a+] = (5jj, i,j = l,---,r (3.1) 

and hence generate the Lie algebra W{r) made out of r copies Wi of the Heisenberg-Weyl 
algebra: 

r 

W{r) = ^W^, (3.2) 

i=l 

where Wi is the linear span of {ai,af , 1}. Denote by |ni, • • • , Uj.) = \ni) ■ ■ ■ \nj.) the normalized 
simultaneous eigenvectors of the r number operators Ni = afai: 

afai\ni,---,ni,---,nr) = ni\ni,---,ni,---,nr), n^ G N, i = l,---,r, (3.3) 

(mi,--- ,mr\ni,--- ,nr) = 5mi,ni ■■■Smr,nr- (3-4) 

The algebra W{r) is represented in the standard way in this number state basis: 

Oj |ni, - - - ,ni,--- , Ur) = y/ui Im, - - - , nj - 1, - - - , n^) , 

o+ |ni,--- ,ni,--- ,nr.) = VnTTTlni,--- ,nj + 1,--- ,nr) . 

The Bargmann realization in terms of coordinates Zj, i = 1, - - - , r in C has 

Q z"i - - - z'*'' 

ai = -^— , a+ = Zj, i = l,---,r, and (zi, - - - ,2:r-|m, • • • ,"-r) = / , ^=- (3.6) 
0-2j Vni!---nr! 



(3.5) 
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4 The Model 

Consider the set of r non-Herniitian operators Hi, i = 1, . • • • , r, defined as follows in hl{W{r)), 
the universal enveloping algebra of W{r): 



H 



/3,c 



O^flfc 



".^E^Mr^^?;(i 



fe=i 



•^ rt (^ ~ ^i) 



J + 



^a+afc + /3 , i= !,■ 



\fc=i 



(4.1) 



Let us observe, especially using the holomorphic representation of Bargmann, that the opera- 
tors X^fc=i a^afcjOi) ^i^) ^i*" X^fc=i Ofc Ofc)^ = !)■■■ j'^j making up the Hi, generate a Lie algebra 
isomorphic to that of the diffeomorphisms in C spanned by (generalized) vector fields of the 
form Z = Yll=ifi{^)7r' + di^ with z = {zi,--- ,Zr). In the coordinate realization where 
Oj = -y=(xi + Q^) and af = -/^{xi — q^), i = 1,- ■ ■ ,r, the operators Hi are third order differ- 
ential operators. However, they can be considered as Hamiltonians, as we shall see below. Now 
the key observation is that they can be simultaneously diagonalized by the multiple Meixner 
polynomials. 

Indeed, consider the states |x,/3, c )) defined through the following combination of the states 
Im,--- ,nr): 



X 



(3,c))=Nlf,^,Yl 



./3,c 



\/nil ■ ■ -rirl 



Im,--- ,nr) , 



(4.2) 



where A''^ o^are constants. Acting with Hi on \x,f3)), we find (with some indices suppressed): 



+ 



n=0 
r 

El 

.k=l 



Uk 



+ 



Ck I- Ci 



/J+E 



nk 



k=l 



|ni,--- ,nr) 



+ /3 + Vrifc V- — ^—-^y^n~TT\ni,--- ,ni + l,--- ,nr)\ 



kT. 



1 



n=0 

r 



1 ^/nll ■ ■ -Hrl 



Mri+gSx) + 



^(^ + N) + Ei 



nk 



k=l 



Ck 



MJx) 



'^yZn — ^—T2'^j(f^ + \"'\ ~ l)^n-e,(2;) \ |rii, • • • ,nr) 

^ (1 - cjY J 



Hence, appealing to the recurrence relations ()2.3p . we have 

Hi'''\x,(3,c)) =x\x,f3,c)) 



(4.3) 



(4.4) 



/3,c 



for all i = 1, • • • ,r. We thus find that in spite of being non-Hermitian, the operators H^' all 
have the same real spectrum. (An explanation is given below.) 

8 c 

Somewhat surprisingly, the operators H^' do not commute. In effect, one finds that 



^^,C jjl5,C 



Qi — a,- + 



(1-q)(1-c,) 



/3 + E« 



fcfflfe 



(4.5) 



fc=i 
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How can this be reconciled with the fact that they can be simultaneously diagonalized? In this 
respect, one observes that 





■^/3,c-^|/ 


' k,/3,c)) 


n=0 


+ (^_^^)(/_^^.)(/3 + H)V^K,...,n.)} 


CXD p 


n=0 


M.,,(.)-M,,., + ^^_;^^^^_^^^ 


to find that 




'H^rc^H^rc 


|^,/3,c))=0 



, lir 



^/nj\nir-- ,nj -I,--- ,nr 



{p + \n\)M^{x) 



Pi,--- j^-r) 



(4.6) 



(4.7) 



from (j2.4p . This shows that the operators H^'^ commute when acting on the states \x,f3,c)). 
These operators therefore define a "weakly" integrable model. 



5 Conjugation operators and isospectrality 

First observe, using (|2.4p in (j2.3p . that M^(x) satisfies "non nearest-neighbor" recurrence rela- 
tions of the form 



xMn{x) = M.n+sSx) 

r 



c, ^^ l« 



1 - C; 



1 -a 



MJx 



+Et 



Cfc 



k=i 



Cfc 



-".fcM, 



n+e — efc 



(x) 



Cfc 



(5.1) 



-nfc(/3+ |n| - 1)M^_^^, i = 1, ■ 



fc=1 (1-Ci)(l-Cfc) 
Consider another set of r Hamiltonians H^''^ defined as follows: 



H 



at + 



l-o. 



■/? + 



k=l H «fc 



-Er 



fc=i 



-Sti 



Cfc 



^ (1 -Ci)(l -Cfc) 



1 - Cj 



Cfc + 



Cfc 



(5.2) 



i = 1,- • • ,r. 



fc=i 



In a manner similar to the procedure followed in Section HI one can easily show that the operators 
H^''^,i = 1, • • • ,r are also diagonalized with the help of the multiple Meixner polynomials. In 
fact one finds that they share with H^''^ the same eigenstates: 

H^'''\x,l3,c)) =x|x,/3,c )). 



(5.3) 



Here it should be noted that -fff '^ 7^ -f^i except for r = 1 , although they have common spectra 
by construction. 

Quite remarkably, each of the operators H^'^ can be obtained from the standard harmonic 
oscillator Hamiltonians in r dimensions 



^0 = Z] "/'^J 



(5.4) 
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by a conjugation. Let Si, i = 1,- ■ ■ , r be the operators defined by 

Cfc 



Si = exp [(1 - Ci)ai] JJ exp 



(1 -Ci)(l -Ck) 



at{f3 + Ho) 



(5.5) 



Using the Baker-Campbell-Hausdorff formula 

exp[X] Y exp[-X] =Y+[X,Y] + ^ [X, [X, Y]] + ■ ■ ■ (5.6) 

and the fact that 

[a+(/3 + Ho), Ho] =at{l3 + Ho), (5.7a) 

[a+(/3 + Ho),a+{/3 + i^o)] = 0, (5.7b) 

it is seen that 

Hi''' = SiHoS-' (5.8) 

with the proper choice for /3 and c in Si implicit. This explains why each of the operators H^'^ 
and hence each of the H^''^ has a real spectrum identical to that of Hq. 

6 The special case r = 1 and SU{1, 1) 

When r = 1, the multiple Meixner polynomials of the first kind reduce to the standard monic 
Meixner polynomials. It is well known that these are related to the discrete series representation 
of the Lie algebra SU{1, 1) which is bounded from below [U EJ [lO]. It is of interest to make the 
connection with our model. 

Let Jo and J± denote the generators of SU{1, 1) satisfying the commutation relations: 

[Jo, J±] = ±J±, [J+, J-] = -2 Jo. (6.1) 

The Casimir operator C is given by 

C = Ji-^{J+J- + J-.U). (6.2) 

We shall take it to be of the following form 

C = d{d-1) (6.3) 

in irreducible representations. The discrete series model of SU{1, 1) which is bounded from 
below is spanned by basis vectors {\d,n) ,n = 0, 1, • • • } on which the generators act according 

to [n] 

Jo \d, n) = {d + n) Jo \d, n) , (6.4a) 



J± \d, n) = y/{d + n± l){d + n) - d{d - 1) \d, nil). (6.4b) 

Consider now the following "Hamiltonians" H in the Lie algebra of SU{1, 1): 

H = J^ + (^) ^0 + .^^ .2 -^+, < c < L (6.5) 
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Using the Baker-Canipbell-Hausdorff formula (]5.6p . it is readily seen that 

H = SJqS-^ (6.6) 

with 



S = exp[(l — c)J-] exp 



' :J. 



(6.7) 



;i-c 

Hence H has the same spectrum as Jq. In order to diagonalize H, one introduces the vectors 

^ y/nl{2d)n 
The eigenvalue equation 

H\d,x) = {d + x)\d,x) (6.9) 

is then readily seen to require that the Pn{x) satisfy the 3-term recurrence relation which is 
obtained from (j2.3p when r = 1, as long as 



(3 = 2d. (6.10) 

This establishes the connection between 5C/(1, 1) and the standard Meixner polynomials. 

Obtaining harmonic oscillator models for these polynomials therefore amounts to providing 
embeddings of SU{1, 1) in the universal enveloping algebra of VF(1). One such well known 
realization corresponds to the metaplectic representation of SU(1, 1). One takes 

J- = y, J+ = -^^ Jo = -(aa+ + a+a). (6.11) 

It is straightforward to check that the defining relations (|6.ip are satisfied, given the canonical 
commutation rule [a, a"*"] = 1. In this realization, the Casimir operator C takes the value 

C = d{d-l) = -^, (6.12) 

lb 

implying that only two series, respectively corresponding to d = j and d = ^ are thus con- 
structed. In other words, this only provides an algebraic model for the standard Meixner poly- 
nomials Mn with /3 = 2 or /? = |. An extension to more general value of /3 can be obtained [T2] 
by letting a and a+ in (j6.1ip stand for annihilation and creation operators of parabosonic oscilla- 
tors, that is by considering a deformation of the Heisenberg algebra defined by the commutation 
relation 

[a, a+] = 1 + 2^ii (6.13) 

where R is an involution operator satisfying B? = 1. In doing so, we construct the SU{1, 1) 
discrete series with d{= ^) equal to ^ -|- | and to ^ -|- | and models for Meixner polynomials 
with the corresponding /3. In any case, we have found that the scheme corresponding to (j6.1ip 
does not extend to r dimensions or to r oscillators, so as to provide an interpretation of the 
multiple Meixner polynomials. 

Now, if we specialize the model presented in section UJ to the case r = 1 (and suppress the 
sole index 1), we first find that the operators H and H coincide and read 

H = H = a+'^-^^ [a+a + + ^^^ [a+(/3 + a+a)] - ^. (6.14) 
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Moreover the operator that conjugates Hq = a~^a into H = H 



IS 



S = exp[(l — c)a] exp 
Hence if we take 



c 



(1 - c)^ 



-a^{l3 + a~^a) 



(6.15) 



J- = a, J-^- = a'^ {(3 + a'^ a) , Jo = a^a + —, (6.16) 

we see that i? + ^ is in the form ()6.5p and that S is the same group element as in (j6.7p . It 
is then straightforward to check that the operators defined in ()6.14p indeed verify the SU{1, 1) 
algebra commutation relations. Furthermore, the Casimir operator C is readily calculated to be 

confirming that d = f and that this embedding of SU{1, 1) in U{W{1)) provides an algebraic 
model for the Meixner polynomials with arbitrary values of /3 > 0. To our knowledge, this 
interpretation of the standard Meixner polynomials through the oscillator realization ()6.14p of 
SU{1, 1) is new. The derivations of properties of the multiple Meixner polynomials stemming 
from our model would hence also be original in the case r = 1. (See the next three sections.) Of 
course, for r > 1, the multi-oscillator generalization of the operators ()6.14p leads to generators 
that realize an infinite Lie algebra isomorphic to Diff (C) as we have already pointed out. 

7 Ladder operators and step relations 

In the remaining sections, we shall illustrate how the algebraic model, presented in section [U 
provides a useful framework to derive properties of the multiple Meixner polynomials. We will 
first focus on the step relations and therefore start with an examination of ladder operators. 
Introduce 

Xi = ai + -^^{/3 + Ho), i = l,---,r (7.1) 

1 - Cj 

and 

k=i ^ 

It is straightforward to check that these operators satisfy the following relations: 

Hr'Xj-X,H^ = -X, (7.3) 

and 

H^-^Y - YH^ = Y, (7.4) 

for i, J = 1, • • • , r. From ()7.3p and ()7.4p . it follows that 

X^\x,f3))=C^,x,^\x-l,|3 + l)), (7.5) 

Y\x,^))=r,.,p\x + 1,^-1)), (7.6) 



where Ci^x,i3 and r/^^/j are proportionality factors to be determined. Therefore Xi (resp. Y) maps 



the eigenstate of -ffj^ with eigenvalue X to the eigenstate of ffj^ (resp. H^ ) with eigenvalue j;+ 
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1 (resp. X — 1). It should be stressed that these operators (much Uke supersymmetry generators 
in quantum mechanics) relate the eigenstates of different operators. Using the definition of 
\x,f3)) and a procedure by now familiar, it is found that (|7.5p translates into 



n: 



mLAx) + T^^(/3 + |n|)M|(x) 



n+ei^ 



1 -C; 



/3+1/ 



ei,.,/3iv;_i,^+iM„r^(x-i). 



(7.7) 



Given the monicity of the polynomials, i.e. that the coefficient of the leading term is one, we 
must have 



^i,x,/5 



N. 



x,P 






[X + a), 



where a is a constant. This gives 



r/3 



/9+1/ 



M^^,^ix) + -^(/3 + \n\)M^{x) = {x + a)M^+\x - 1) 

i Ci 



(7.8) 



(7.9) 



In order to determine the constant a, we set n = 0. We immediately find that a = since we 
know from the recurrence relations (12.31) that 



mS^'^ = X 



-/3. 



1- Ci 

We thus arrive at the following backward relation 

xM^+\x - 1) = m|^. (x) + -^^(/3 + |n|)M|(x), i 

Similarly, it is seen that ()7.6p amounts to 



, r. 



K,p 



k=i " 

The monicity of the polynomials here implies that 



V.,^K+i,p_,Mt\x + l) 






thereby giving the following forward relation: 



(7.10) 



(7.11) 



(7.12) 



(7.13) 



M, 



t\x + 1) = M|(x) + J: r-^Ml,^ix 

11^ ^k 



k=l 



(7.14) 



8 Intertwiners and contiguity relations 

Contiguity relations can be obtained in a fashion similar to the way the step relations were 
derived in the last section. One observes that the operators 



X^ = a, + - W + Ho), i = l,---,r 

1 - a 



(8.1) 



and 



Yi 



'-^T 



Ck 



k=l 



Ck 
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satisfy the following intertwining relations: 



and 



HrX, = X,H^ 



H^-^Y = YH^, 



for i, j = 1, • • • , r. From ()8.3p and ()8.4p . we conclude that 

Y\x,P))=fi,,)3\x,P-l)), 



(8.3) 



(8.4) 



(8.5) 
(8.6) 



with ^i^x,/3 s-nd 'qx,p proportionality factors again. Note that in contradistinction to (j7.5p and 
(j7.6p . the eigenvalue x is unscathed in (j8.5p and (j8.6p with Xi (resp. F) mapping the eigenstate 
of H^ with eigenvalue x to the eigenstate of i?f (resp. -fff ~ ) with the same eigenvalue x. 
We now proceed as in the last section. Eq. (j8.5p is seen to imply that 



K,p 






-(/3 + |n|)M|(x) 



/3+1/ 



?i,.,/3iv:_i,^+iM„r^(x). 



(8.7) 



with monicity forcing 



SJ.a;,, 



iv: 



iv: 



-{x + d), 



a;,/3+l 



where a is a constant. Here also a is determined by taking n = and one finds a = f3. One thus 
arrives at the following relation 



{x + /3)m|+^(x) = m|^..(3^) + r. (/3 + |n|)M|(x), i = 1, • • • , r. 

A second relation is obtained from ()8.6p which entails 

r 

Cfc 



(8.9) 



1 1 f^ 



(8.10) 



k=l 



with 77^.^^ required to be 



Nl 



Vx,l3 



x,P 



N 



.11) 



x,l3-l 



Here it should be mentioned that combining ()7.1ip and (j8.9p . one finds 

xM§+\x - 1) - q(x + /3)m|+1(x) = (1 - c,)m|^. (x), 



for i = 1, • • • , r and using (j7.14p and (|8.10p one gets 

r 

m|(x + 1) - m|(x) = ^nfcM|+4(x). 



.12) 



.13) 



ifc=i 
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The raising and lowering relations (|8.12p and (J8.13P have been shown [13] to yield a difference 
equation for the multiple Meixner polynomials that reads 

r 

n [cj{x + /3) - xT-i] o (T, - I)M^/{x) 

(8.14) 
= ^(cfc - l)nk n [ci(x + /3) - xT-^] M^/{x), 

fc=l j^k 

where / is the identity operator and T^ is the shift operator for x, i.e. Tx[f{x)\ = f{x + 1). 

We can also derive another "difference" equation for the multiple Meixner polynomials not 
in X but in (3 by combining the relations (18. 9p and (18.100 : 



n [{cj - l)(x + /3 + l)Tp + (/? + \n\)] o {Tp - I)M^ 



^^''x) 



^ ^ (8.15) 



r 



= Y. ^'^^'^ n [(^J - 1)(2; + /3 + 1)7> + (/3 + |n|)] m|+^''^(x), 
fc=i jyfc 

with T^ the shift operator for 13. 

9 Generating function 

A generating function for the multiple Meixner polynomials of the first kind is readily obtained 
from the model with the help of the Bargmann realization. With Oj and af realized as in ()3.6p . 
the Hamiltonian H^''^ reads 

H^^'=j^hu{z)^ + g^{z). (9.1) 



k=l ^'^ 



where 



fik{z) = 5jfc + ( -r^. + ^_^ ) ^fe + X] (1 _^c.')2 ^j^fc' (^-^^ 



1 - Ci 



The eigenvalue equation -fff ''^ I x, /?)) = x | x, /?)) thus becomes a first order linear partial differ- 
ential equation in C. One can check directly that it has for solution the function 

with the constant of integration chosen so that ^xi^^ = 1- Recalling how the number states 
are represented, see (13. 6p . in the Bargmann realization we have: 



(zi,---,z.|x,/3,c)) = iV;^^ ,"..V. ^"'---^"^- (9-^) 



ni! • • • n^.! 

n=0 
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Note that (0, ••• ,0|a:,/3,c)) = N^ a] hence, (fx{z) = jj^^ {zi,- ■ ■ , Zj,|x,/3,c )). Bringing this 
together gives the generating function formula 




10 Conclusions 

Let us summarize our results. As was done previously for the multiple Charlier polynomials, 
we have provided here an algebraic interpretation of the multiple Meixner polynomials of the 
first kind. They were shown to arise in the construction of the common eigenfunctions of a set 
of r non-Hermitian Hamiltonians H^''^ i = 1, • • • , r. These turn out to be in involution only 
"on shell" where they define a weakly integrable system. This model is not a multi-variable 
extension of the familiar connection between the one-dimensional harmonic oscillator and the 
Meixner polynomials via the metaplectic representation of SU{1,1). The isospectrality of the 
Hi has been traced back to the fact that their spectra are identical to those of operators Hi 
that can be obtained from the standard r-dimensional harmonic oscillator Hamiltonian through 
conjugation by some specific operators. Finally, the model has been shown to give interpretations 
to the step and contiguity relations as well as to a generating function. 
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